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We study the phenomenology of fixed-target elastic J/tp photoproduction in the NRQCD factor- 
' ization formalism. Our the goal is to test an essential feature of this formalism — the color-octet 

. mechanism. We obtain an order-of-magnitude estimate for a certain linear combination of NRQCD 

' color-octet matrix elements. Our estimate is consistent with other empirical determinations and 

T— I , with the u-scaling rules of NRQCD. 
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O ■ I- INTRODUCTION 

We study the phenomenology of the NRQCD factorization formalism, which is a framework for writing inclusive 
I rates of production and decay for quarkonium. Specifically, we consider photoproduction of J/ip- This process probes 
. an essential feature of the NRQCD formalism — the octet-mechanism. We approach our subject with the question: 
' can existing photoproduction data be used to test the octet-mechanism? 
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(3 ■ II. COLOR-SINGLET APPROACHES TO PHOTOPRODUCTION OF J/i/j 

^ ' 

Before discussing the NRQCD factorization formalism, we first discuss other theoretical frameworks for the calcula- 
tion of charmonium production rates. One of the earliest computations of J/ip production was carried out by Berger 
Qh[ and Jones in 1981 0]. There, the authors give an expression for the rate of production of J/ip in 7-nucleon collisions, 
ji^' as calculated in the so-called "color-singlet" model. In this model, one adopts a nonrelativistic boundstate picture 
I to describe the J/ip calculations of J/ip production are based upon the amplitude for the generation of a cc pair 
. in a color-singlet ^5*1 configuration with small relative momentum {i.e. |q| ^ rric). The cc subprocess considered by 
Berger and Jones is shown in Fig. 1. In it, a gluon from the nucleon "fuses" with the photon to form a hadronizing 
cc pair that recoils against a gluon jet. These diagrams represent the leading color-singlet contribution to 



X 



H ■ + N ^ J/iP + X . (1) 

The color-singlet model formulation of J/ip photoproduction has been studied extensively . These analyses hinge 
on the experimental kinematic parameter 



Pj/4, ■ Pn 
P.-Pn ' 



(2) 



where Pj/ijj is the J/ip four- momentum, is the four- momentum of the initial state nucleon, and is the initial 
state photon four- momentum. (Note that in the lab frame of fixed target experiments, one has z = E^/E^.) It 
is found that once next-to-leading order QCD corrections to the Berger- Jones results are taken into account, the 
color-singlet model can adequately explain the experimental data for the kinematic regime pT > IGeV and z < 0.8, 
where px is the transverse momentum of the J/ip in the center-of-mass frame. 
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FIG. 1. Leading order diagrams for the photoproduction of a cc in a color-singlet, ^5*1 state. 



However, the color-singlet model alone cannot explain the total J/iIj photoproduction cross section. It severely 
underestimates the rate of production in the region z > 0.9 Q [H, which (irnportantly) is the region containing nearly 
all of the events. In production of J/V" in this high-z range, the "X" in Eq. |l| represents either an elastically scattered 
nucleon, or, a nucleon (or resonance thereof) and light hadrons that have very low energy in the rest frame of the 
final state nucleon. Such events, sometimes called "diffractive scattering," have been modeled in Refs. |^ and [Q. 
The essential features of the diffractive mechanism are the creation of a heavy quark pair from the incoming photon 
and the exchange of some color-singlet gluon combination between the cc pair and the proton. (Diffractive processes 
can in general entail the creation of a heavy quark pair from an incoming gluon also.) In the analysis of Refs. ||^ and 
, the amplitude for the completely exclusive process 7p J/ipP is factored into three pieces: the process 7 — > cc, 
the scattering of the cc system on the proton via (colorless) two (or multiple) gluon exchange, and the formation of 
the J/ip from the outgoing cc pair. The exchanged gluons are taken to have very low transverse momentum. This 
mechanism is considered in Ref. to serve as a probe of the gluon density in the proton, the cross-section being 
proportional to [xfg^pf^x)]'^: 

da{^p J/il:p)/dt cx [xfg/Nix)]^ (3) 

where a; is a dimensionless ratio of kinematic variables. Indeed there is remarkably good agreement between the 
formula of Refs. |6| and [Q and high-z fixed target and HERA data (see Fig. 2 of Ref. 0]). 

The photoproduction cross section can be described at low z by the color-singlet model, and at high z by the 
approach of Refs. and 0. It is interesting however to consider the possibility that the total J/ip photoproduction 
cross section can be parametrized within a single coherent theoretical framework. Such a framework is provided 
by the nonrelativistic QCD (NRQCD) factorization formalism of Bodwin, Braaten, and Lepage This formalism 
is intended only for the parameterization of inclusive production of quarkonium. Here, "inclusive" means that a 
perturbative partonic process produces the heavy quark pair and (possibly) some other specified final state partons; 
the heavy quark pair hadronizes into a specific quarkonium boundstate but no mention is made in the formalism of 
the fate of the other partons regarding the baryons into which they hadronize; all possibilities are included] moreover, 
a large number of light quanta (generally pions) are expected to accompany this inclusive production. 

At the outset, the choice between the diffractive scattering formalism and the NRQCD factorization formalism 
hinges on the questions that one asks, that is whether one desires to model exclusive or inclusive events. We will see 
below however that the line between these two cases is blurred by the fact that in the set of all data, a large portion 
of the J /ip photoproduction events have a certain exclusive character — that is, the final state contains very few of 
the extra quanta (such as pions) that one would expect in inclusive production. 

III. NRQCD FACTORIZATION FORMALISM 

According to the NRQCD factorization formalism, the inclusive cross section for quarkonium production is expressed 
as a sum of products known as a "factorization formula." For the photoproduction of J/'0; the factorization formula 
is 
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a(7 + ^ - J/V^ + ^) = E -^m'/^nM . (4) 

_ Tflr 
n ^ 

In the above expression, the index n labels the initial color and angular-momentum quantum numbers of the cc pair 
produced by short-distance physics. The short-distance coefficients F{n) contain only effects of distance scales of 
order 1/toc (where rric is the charm quark mass) or smaller; they can be calculated, using Feynman diagrams, as 
a perturbative expansion in asirric)- Effects of longer distances, including effects related to the hadronization of 
the cc pair into the J/ip boundstate, are parametrized by the NRQCD matrix elements (0|C'''/'''(7i)|0); their relative 
importance can be determined using the NRQCD u-scaling rules where v is the typical relative velocity of the c 
and c in the J/t/j boundstate. c?„ is an integer associated with the energy dimension of the operator O'^^^in). 

The NRQCD formalism has several features that make it highly compelling. Firstly, it is rigorous in that it is 
based upon an effective Lagrangian for heavy quarks derived directly from full QCD. Secondly, it goes beyond the 
color-singlet model in that it takes into consideration processes in which a heavy quark pair is produced initially in a 
color-octet state and attains color-neutrality and hadronizes via the emission and absorption of soft gluons; the rates 
due to such color-octet channels are rigorously and systematically parametrized in terms of NRQCD matrix elements. 
Finally, the NRQCD factorization formalism provides a solution to the problem of the infrared divergences that arise 
in naive color-singlet-model calculations of P-wave quarkonia decay and production pO| | ; in the NRQCD framework, 
these divergences are absorbed into appropriate (color-octet) matrix elements. 

The factorization hypothesis expressed in Eq. ^ is considered to be most valid when applied to differential cross 



sections with pT much greater then the hadronic scale 1 GeV |29|. However, it is not yet clearly understood in the 
NRQCD theoretical community whether total cross sections can be expressed in the factorized form. In this paper 
we will assume that they can be. If our assumption is valid then one expects relative corrections from higher-twist 
operators to Eq. || to be 

O ((1 GeYf/ml) , (5) 

with rrir in GeV. 



IV. NRQCD PREDICTION FOR PHOTOPRODUCTION 

Let us now turn to a discussion of the NRQCD prediction for the total cross section for photoproduction oi J/^. To 
write down the leading pieces of the NRQCD prediction, we must first determine which are the most important terms 
in the factorization formula. The numerical sizes of the NRQCD matrix elements can be estimated by determining how 
they scale with v, which, as stated earlier, is the typical velocity of the heavy quarks in the quarkonium boundstate. 
(w^ ~ 0.25.) Combining u-scaling estimates of {0\O'-'^'^ {n)\0) with the ag-scaling of the F{n), it is possible to determine 
the relative importance of the terms in the factorization formula, Eq. ^, in regard to the double expansion in v and as . 
One finds that the leading contributions to the total photoproduction cross section come from the production of a cc 
pair in a color-singlet ^Si state (the Feynman diagrams are shown in Fig. 1), and from the production of a cc pair in a 
color-octet ^^o, "^Pq, or '^^2 state (the Feynman diagrams are shown in Fig. 2). The leading color-singlet contribution 
to the rate is proportional to a^v^, while the leading color-octet contributions to the rate are proportional to a^v"^ . 
Since as{mc) , we expect the two effects to contribute roughly with the same size to the total cross-section. 

Note that in the case of J j'^ production where the cc pair is produced initially in a color-singlet state (as in 
Fig. 1), there must necessarily be a final state hard gluon, for color conservation. If this hard gluon is of relatively 
low (high) energy, then the J/-)/) is produced in the high-z (low-z) region. For the color-singlet contribution, we thus 
see that there exists some continuous z-distribution. On the other hand, as to the case (in Fig. 2) in which the cc 
pair is initially produced in a color-octet ^5*0, ^Pq-, or ^P2 state, there is no hard final state gluon at leading order; the 
color-octet mechanism contributes therefore only to the high-z region at leading order. The naive picture is that the 
z-distribution from this sort of contribution is an idealized Dirac function at z = 1. We will refer to the color-octet 
process in Fig. 2 as "forward-octet." 

We now present the color-singlet J/V' photoproduction rate, the underlying process of which is shown in Fig. 1. The 
NRQCD result for the color-singlet contribution is simply proportional to the color-singlet model expression |Q] |^ 

rfacs ^ 647regaemas(^)^m,|j;,(0)p s^(s - 4mg)^ + e{t - 4mg)^ +u^{u- Amlf 
dt ~ 3s2 (s-4m2)2(t_4m2)2(u-4m2)2 ^ ' 

where as{t) is the strong coupling constant evaluated at the typical scale of the interaction, i, where Uem is the 
electromagnetic coupling constant, where is the fractional charm quark charge, and where s, t, and u are the usual 
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FIG. 2. Leading order diagrams for the photoproduction of a cc in a color-octet state with angular momentum configuration 
^5*0, ^Po, and ^P2. These dia grams are the underlying subprocess for the forward octet contributions to J /ip photoproduction. 



Lorentz invariant Mandelstam variables. Here Rs{0) is the radial color-singlet wavefunction evaluated at x = 0. It is 
related to the color-singlet ^ Si production matrix element through 

(0|<^f5i)|0) = |-|i?,(0)p(l + O(^;2)) . (7) 

We next present the leading color-octet ("forward octet") photoproduction rate, the underlying process of which is 
shown in Fig. 2. The factorization formula for the forward-octet contribution can be written as follows: 

mi 

TO* trq. 

One can lift (with the help of a suitable color-factor replacement), the short-distance coefficients F(8, ^iSq) and 
F(8, ^Pj)^ from the results for the process gg J/ip given in Ref. One obtains 



^^(8, ^5o) = 5{^ml s) 

nrir 



F{8, ^Po) = 37r^ega.(2m.)ae,n _ 
rric 

F8(8, 3P.) = l ^^^c«s(2m.)a..^ 

5 TOc 

Inserting these into the factorization formula in Eq. ^, we obtain the subprocess cross-section 

o'(75 J/^ + ^)fo = — 5{Ami - s) 6 , (10) 



where is given by 



o ^ ml'\'s,m + J-^{Q\oi'\^p,m . (n) 

mi 



The above expression is derived with the help of the relation 
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{Qpi'^i^pjm = (2j + imoi'^i^p^m [i + o{v^)^ . (12) 

We next convolute the subprocess cross sections given in Eq. ^ (color-singlet) and Eq. ^ (forward-octet) with the 
gluon distribution function to obtain the leading NRQCD factorization formalism prediction for the total photopro- 
duction cross section: 

a{-iN J/V' + X) 

= f dx fg/N{x) (crcsiig J/ip + X) + apoiig J/'(p + xj^ 

dx fg/N{x) 7ras(2TOc) B S (Ami - + 



2 2 
ml 

* 27s2 (s -4m2)2(t-4m2)2(u-4m2)2 ]' ^ ' 

where s, t, and u are the usual Mandelstam variables for the subprocess cross section, and fg/N{x) is the gluon 
structure function. 



V. COMPARISON TO EXPERIMENT: CAVEATS 

The leading NRQCD factorization formula for the total photoproduction cross section (Eq. |l^) depends on the two 
phenomenological quantities: (0|of '^(^S'^jO), and 9 = {Ops^''' So)\0) + :^{0\O^^''' Po)\0) . The first parameter 
(the color-singlet matrix element), is well determined from measurements of the decay rate of J/ip to two leptons [p^ . 
As to the second parameter (the forward-octet combination O), the matrix elements contained therein are poorly 
constrained thus far, and so a testable prediction of the photoproduction cross section is not yet possible. 

Given this state of affairs, we propose to use the factorization formula in Eq. ^ to obtain an estimate of the value 
of Q via a comparison with fixed target photoproduction data. We find that a determination of Q is complicated by 
the manner in which J/tp photoproduction data is presented in the experimental literature. Typically, measurements 
of the total inclusive cross section are not made; thus, our Eq. |l^ — a total inclusive cross section — cannot be 
associated in an unqualified straightforward manner with some published empirical result. We must contend with the 
fact that data is taken in either the inelastic or elastic regime. The inelastic regime is generally considered to be the 
region where z is below 0.8 or 0.9. On the other hand, the conventional definition of the elastic regime is not at all 
well established; it is in general considered to be the region near z — 1, with the data consisting of those events in 
which is detected a muon pair with invariant mass along with hadrons having up to about 5 GeV of additional 
energy. To test the forward-octet feature of NRQCD (which is our goal), it would seem at the outset that a good 
strategy would be to select some range of z (from a lower bound around 0.8 up to unity) for which to calculate the 
total (color-singlet plus color-octet) theoretical prediction, and to compare that result to the elastic data. However 
the contribution of the color-singlet term is quite negligible for high z, and is actually found to be nearly an order 
of magnitude below the data Q ||]. The unimportance of the color-singlet contribution dacs/dz at high z and the 
vagueness surrounding the correct lower bound of z to use when integrating over z in order to compare to elastic data 
suggest that we might as well neglect the color-singlet contribution altogether. Why not simply go ahead and compare 
the elastic data to the forward-octet contribution? Before we proceed, some reservations must first be discussed. 

A major reservation concerns the issue of "z-smearing." Although we ultimately do not take this stance, it is 
arguable that one actually does not expect the forward-octet piece to contribute significantly to photoproduction 
events for which z is nearly unity. Although the subprocesses described in Fig. 2 appear, at a first glance, to create 
J ftp particles with total energy close to that of the incoming photon {i.e. with z — {Ej/^,)/{E^) = 1), this is not 
exactly the case since the cc system must emit or absorb gluons in order to make the transition to a color-singlet state 
and hadronize into a J/ip boundstate. Such gluons typically have energy and momenta of order mcW^ in the rest-frame 
of the cc system . The presence of these gluons implies that the ratio of the energy of the final charmonium particle 
to that of the initially produced heavy quark pair is not exactly unity but, rather, is expected to be < z 1 — 0{v^). 
The z-distribution is smeared away from the idealization of the Dirac delta S{1 — z). Similar reasoning has been given 
in In fact, even in quarkonium production events involving color-singlet channels, one still might expect such 
a disparity between the energy of the initially produced heavy quark color-singlet pair and the final boundstate; the 
heavy quarks emit pions and other quanta so as to adjust their energy-momentum to that of heavy quarks in the 
boundstate. 
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The above arguments concerning z-smearing seem reasonable, but experimental information challenges them 
severely. In fact, quite contrary to what one might expect, it turns out that the great bulk of photoproduction 
events actually lie in bins above z = 0.95. For example, of the approximately 850 J jil) photoproduction events col- 
lected by the HI experiment, 700 lie in the bin above z — 0.95 This includes both the events in which the proton 
is reconstructed, and those in which the proton is dissociated. So, assuming that these events are due to the NRQCD 
forward octet mechanism, wc must conclude that the smearing of the z-distribution is not a big effect. |^ 

There is yet another reason for which, at the outset, one might have reservations concerning the fitting of Eq. |l^ 
to high-z fixed target elastic data: the question of the inherently inclusive nature of the NRQCD formalism versus 
the possible lack of complete inclusiveness in the data. However, as it happens, the fixed target data do appear to 
be sufficiently inclusive for our purposes. This is because the HI data as was mentioned just above, show that 
most photoproduction events are associated with z > 0.95. If we can count on all high energy {^/s ^ M^) J/ip 
photoproduction to have this feature, then the elastic fixed target data that we are using should surely include the 
great majority of events, and can be considered inclusive. 

The above arguments support the decision to simply associate the high-z fixed target data with the NRQCD 
forward-octet contribution. 



VI. COMPARISON TO EXPERIMENT: DETERMINATION OF 9 



We now proceed to use experimental data from the E687, NA14, E401, NMC and E516 experiments to evaluate Q. 



We will then show that the value we determine in this way is consistent with HERA (HI ||5| and ZEUS |15|) data. 

In Fig. 3 we display our fit of the forward-octet piece in Eq. ^ to elastic cross-section measurements from the fixed 
target Fermilab experiments E687 NA14 0, E401 |l|], NMC H, and E516 In the same figure we also 
show data from the SLAC experiment |2l|| ; we do not however include this data in our fit since it falls in the low 
energy (i/s ~ M^) regime. As indicated in Fig. 2, we have performed our analysis using CTEQ MRS |^^, and 
GRV structure functions. Wc have taken rric — 1.5 GeV and as{2mc) — 0.26. With this choice of parameters we 
find 

e = 0.02 GeV^ . (14) 

Note that there are large theoretical errors associated with this determination of 8. These will be discussed in the 
next section. 

Having estimated Q from fixed target Fermilab data, it is interesting to see whether we can use this result to make 
a successful theoretical prediction for the higher energy elastic data from HERA. The prediction is shown in Fig. 4, 
where it is seen that our extrapolated theoretical curve (generated with O — 0.02GeV'^) is indeed consistent with 
HERA data. On this plot we show data from the experiments HI M and ZEUS |15|. (Also shown on this plot are 



measurements from the Fermilab experiments E401 [|l^ and E516 ||20||, which are a class of data that the HERA 
collaboration includes on their own plots for comparison purposes. Note that the E401 and E516 data appearing in 
this figure are not the same as those appearing in Fig. 3.) 



VII. THEORETICAL UNCERTAINTIES 



Here we discuss the sources of theoretical error bearing upon our determination of Q. We have already stated that 
one expects an error do to the neglect of "higher-twist" operators of O ((1 GeV)^/m^). Apart from this limitation, 
there are other important sources of theoretical uncertainty. These are associated with as corrections, with relativistic 

corrections, and with the value of the parameter rric. We now address each of these issues. 

As to as and relativistic corrections, to determine their size, one really must calculate them. When one is not 
in possession of such results, the order-of-magnitude of the corrections can nevertheless be estimated by examining 
similar processes where next-to-leading order corrections have been computed. In this matter we are fortunate, 
because there exist calculations of the ag corrections Q and corrections |p5| to the color-singlet term, Eq. ^ The 



^It must be mentioned in passing that in the exclusive diflractive scattering, as well, one would expect z — smearing to occur. 
The c and c, which emerge from the diffractive scattering in a color-singlet state, would have to readjust their momenta to that 
of quarks in the quarkonium boundstate, and this necessarily entails the emission of soft quanta of momentum of the order of 

2 

mv . 
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FIG. 3. Fit to fixed-target experiments. Shown are data points for tlie elastic photoproduction cross section measured at 
fixed-target experiments. Also shown are the theoretical curves evaluated numerically using five different gluon distribution 
functions. In the curve is adjusted to best fit the data. 
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FIG. 4. Comparison of NRQCD forward octet theoretical prediction to HERA data using tlie value of Q determined from 
fixed target experiments. Shown are data points for the elastic scattering cross section measured at various HERA and Fermilab 
experiments. 
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extrapolations of these results to a color-octet production process must be considered no better than educated guesses; 
we seek only order-of- magnitude guidance. According to Ref. the as corrections to the color-singlet piece increase 
the theoretical result for that part of the rate by roughly a factor of two; we expect Ois corrections of similar size 
therefore for the forward-octet terms. Concerning relativistic corrections on the other hand, the results of Ref. [2^ ] 
show that corrections to the color-singlet part of the rate contribute roughly an additional 50%, and so we expect 
the same for the forward-octet part of the rate. 

As to the theoretical error due to uncertainty regarding the value of the parameter rric, we estimate this by varying 
rric over the range 1.3 GeV to 1.7 GeV. We find that over this range, the cross section varies by a factor of 4. 

Clearly, in light of the large size of the various theoretical uncertainties, the value we determine for 8 can only be 
regarded as an order-of-magnitude estimate. 

The uncertainty due to the gluon structure function is negligible in our analysis. We observe that different standard 
parameterizations of the gluon structure function yield results that vary little compared to the relatively large errors 
discussed above. Therefore, it is not possible to make any conclusions about the gluon content of the nuclcon in the 
present context. 



VIII. COMPARISON TO OTHER WORK 



We now confront other measurements of the NRQCD J /ijj production matrix elements with our measurement of O. 

Beneke and Rothstein ||l^ calculate the leading terms in the NRQCD factorization formula for a{TrN — > J/ip + X). 
This hadroproduction reaction involves — among other processes — the forward octet process in Fig. 2, but with 
two incoming gluons instead of the photon and gluon as in photoproduction. As a consequence, (0|C'g^'''(^5o)|0) and 
(0|Og^'''('^Po)|0) appear in hadroproduction rates in the exact same linear combination (0) as they do in photoproduc- 
tion. Comparing their factorization formula to data, Beneke and Rothstein obtain 6 = 0.03 GeV^. This is consistent 
with our result, given the large theoretical uncertainties. 

A different linear combination of (0|Cg^'''(^5'o)|0) and (0|C'g^'''('^Po)|0) is determined from an analysis of J/tp 
production at the Tevatron . A fit of theory to CDF data for J/ ip produced at moderate pr gives 

{Olo'J^'CSom + ^(0|08^/'^(3Po)|0) = 4.38± 1.15tJ:52 x 10"^ QeV' . (15) 

Here we have quoted the fit done by Beneke and Kramer using the CTEQ4L structure function. Again, within the 
theoretical uncertainty, this result is consistent with ours. 

It is interesting to consider the possibility that our value for Q (determined from photoproduction) is accurate, 
and that the measurement given in Eq. |l5| (determined at CDF) is also accurate. Upon solving Eqs. |lj and |l^, one 

finds that (0|C'g^'''('^Po)|0) is negative. Though this may at first appear to be inconsistent and even meaningless, 
it is actually not unreasonable, if one considers that the matrix elements determined here are renormalized matrix 



elements |27 



The elastic photoproduction calculation carried out in the present paper has also been carried out by Cacciari and 
Kramer Q. The expression they obtain for the forward-octet terms is in agreement with our Eq. |l^. However, in their 
analysis, these authors arrive at conclusions different from ours. Cacciari and Kramer assume that the production 
matrix elements must be positive. They use the result of the analysis of Cho and Leibovich (which is similar 

to the Beneke and Kramer result, given here in Eq. |l^), to estimate the order of magnitude of (0|Og^'^(^S'o)|0) and 
of (0|Cg^'''('^Po)|0). They find that, under these assumptions, the forward-octet term in the theoretical prediction 
for the \ow-pT high-z cross section is one order of magnitude greater than experiment. Their conclusion is that the 
NRQCD factorization formalism might be in confiict with the data. This potential dilemma is alleviated, though, by 
the fact that, as is pointed out above, not all the matrix elements need be positive. 

Cacciari and Kramer also calculate the color-octet contribution to the inelastic (moderate z) photoproduction cross 
section da/dz. (We have not computed the color-octet contribution to this regime of the differential cross-section.) 
Cacciari and Kramer conclude once more that the NRQCD factorization formalism appears at variance with data. In 
this instance, not only is the theoretical prediction several times greater than the data, but, this time, the result of 
the computation of the cross-section (here, the inelastic color-octet cross-section) cannot be affected appreciably by 

the fact that {Op^' Pq)\Q) can be negative . Nevertheless, it is clear from the long list of sources of theoretical 
uncertainty that no strong conclusions can be drawn from this discrepancy regarding the correctness of the NRQCD 
factorization formalism. 
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IX. CONCLUSIONS 



We have calculated the leading order (in coupling constants and v) color-singlet and color-octet contributions to 
the NRQCD factorization formula for the photoproduction of J/tp. At this order the total cross section depends 

on the two phcnomenological parameters: (0|C'i/'^(3S'i)|0), and 9 = (0|Cl8^'^(^S'o)|0) + ^{0\O^^''' {^Po)\0). Though 
the color-singlet matrix element is well determined, the color-octet matrix elements contained in the forward-octet 
combination 9 are poorly constrained. Therefore a testable prediction of the photoproduction cross section cannot 
be generated. 

Be that as it may, it is nonetheless possible to use fixed target elastic photoproduction data to obtain a rough 
estimate of 9. There are some reservations to this approach, though. Firstly, there is the issue of the "z-smearing" 
of the forward-octet contribution. One might be concerned that the forward-octet contribution is centered around 
< z >— 1 — 0{v^), not z — \. Secondly, there is the issue of inclusiveness. NRQCD factorization formulas apply 
to inclusive production, and it would seem that the region z > 0.95 is not inclusive enough. However both of these 
arguments are challenged by experimental data collected by the HI collaboration, which shows that of 850 J/ip 
photoproduction events 700 are contained in the region z > 0.95. 

Proceeding with a fit of the forward-octet contribution to fixed target J/^ elastic photoproduction data we determine 

9 = 0.02 GeV^ (16) 

Given the large theoretical uncertainties associated with this calculation, it is best to regard this determination as an 
order of magnitude estimate. We also show that, using the value of 9 quoted above, we can explain the higher energy 
HERA data. Furthermore the value we determine for 9 is consistent with other determinations of the color-octet 
matrix elements from hadronic J/ip production at fixed target and collider experiments. Moreover, our value of 9 is 
consistent with the w-scaling rules of NRQCD. 

It must be appreciated that our NRQCD forward-octet result for the photoproduction rate in linear is the gluon 
structure function of the proton (see Eq. |l^) while the diffractive rate calculated in Ref. Q and Q is quadratic in the 
gluon structure function (see Eq. |3|). Thus, the two methods are essentially irreconcilable. However, one observes that 
both approaches are capable of convincingly describing data. The extent to which one might "prefer" NRQCD over 
diffractive scattering as the correct explanation hinges on whether one accepts the hypothesis that, when calculating 
the production rate for heavy quarks, to avoid double-counting one ought not include certain additional mechanisms 
that enhance cross-sections such as a diffractive scattering Q . 

We have failed to resolve this issue. Perhaps a valuable diagnostic test to probe the photoproduction mechanism 
would be to compute various polarization rates in both formalisms, for confrontation with experiment. This includes 
polarization of the incoming photon, as well as of the final J/ip. Also, a valuable diagnostic would be to compare 
data to various differential cross-sections calculated in both formalisms. 
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